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ABSTRACT 

Using 10 095 galaxies {B < 20 mag) from the Millennium Galaxy Catalogue, we derive 
B-band luminosity distributions and selected bivariate brightness distributions for the galaxy 
population subdivided by: eyeball morphology; Sersic index (n); 2dFGRS 77-parameter; rest- 
(u — ?■) colour (global and core); MGC continuum shape; half-light radius; (extrapolated) cen- 
tral surface brightness; and inferred stellar mass-to-light ratio. All subdivisions extract highly 
correlated sub-sets of the galaxy population which consistently point towards two overlapping 
distributions: an old, red, inert, predominantly luminous, high central-surface brightness sub- 
set; and a young, blue, star-forming, intermediate surface brightness subset. A clear bimodality 
in the observed distribution is seen in both the rest-(u — r) colour and log(n) di st ribut ions. 
Whilst the former bimodality was well established from SDSS data JStrateva et alT200ll) . we 
show here that the rest-(u — r) colour bimodality becomes more pronounced when using the 
core colour as opposed to global colour. The two populations are extremely well separated in 
the colour-log(n) plane. Using our sample of 3 314 (B < 19 mag) eyeball classified galaxies, 
we show that the bulge-dominated, early-type galaxies populate one peak and the bulge-less, 
late-type galaxies occupy the second. The early- and mid-type spirals sprawl across and be- 
tween the peaks. This constitutes extremely strong evidence that the fundamental way to di- 
vide the luminous galaxy population (A/s^^.^, — 5 log h < ~1Q mag, i.e., dwarfs not included) 
is into bulges (old red, inert, high concentration) and discs (young, blue, star-forming, low 
concentration) and that the galaxy bimodality reflects the two component nature of galaxies 
and not two distinct galaxy classes. We argue that these two-components require two inde- 
pendent formation mechanisms/processes and advocate early bulge formation through initial 
collapse and ongoing disc formation through splashback, infall and merging/accretion. We 
calculate the B-band luminosity-densities and stellar-mass densities within each subdivision 
and estimate that the z « stellar mass content in spheroids, bulges and discs is 35 ± 2 per 
cent, 18 ± 7 and 47 ± 7 per cent respectively. 

Key words: galaxies: fundamental parameters — galaxies: luminosity function, mass func- 
tion — galaxies: statistics — surveys 



1 INTRODUCTION 

Galaxies exhibit remarkable diversity, and as such it is not clear 
how useful collective studies are, i.e., those that bunch all galax- 
ies into a single sample. For instance knowledge of the variation 
of the median galaxy size with redshift has no unique interpreta- 
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tion if multiple evolutionary mechanisms exist. Given the estab- 
lished diversity in luminosity, size, shape, colour, star-formation 
rate, metallicity, gas and dust content, it seems highly likely that 
multiple evolutionary paths, mechanisms and time-scales do exist. 
From a theoretical perspective there are several proposed modes 
of evo l ution: mo n olithic collapse 1 Eggen, Lyndon-Bell & Sanda^ 
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Il986h . satellite accr etion jSearle&ZinnlligVSh. secular evolu- 
tion ( see review bv | K orrnend y^^Kerml£utTl2004l) and splash- 
back jpukueita & Peeble^ ^2003>. The number of secondary (and 
mainly environ mentally depen dent) processes is even higher, e.g., 
tidal formation jS arnes & Hem guist 199 2), ram-pressure stripping 
iGuim & GottI I1972I ). strangulat ion tB alogh. Navarro & Morrij 
GOOOp, harassment ( Moore et alj Il996l). squelching ( TuUv et al. 
I2OO2 1. threshing ( Bekki. Couch & Drinkwater 2001), and cannibal- 
ism ^Dstriker & Hausman 1972), for example. Each may apply, in 
differing degrees at different times with environmental and initial 
condition dependencies. A key question is whether any of the read- 
ily observable properties of the galaxy population today contains a 
clear imprint of these distinct formation mechanisms. One purely 
empirical way to attempt to identify such connections is to look for 
naturally occurring sub-groupings within some region of galaxy pa- 
rameter space. 

At this point the natur al subdivision(s) of g alaxies is not en- 
tirely clear (see for example Isianton et all2003a) and the quest to 
determine a fundamental basis is ong oing and constitu tes a contin- 
uation of the same questions asked bv lHubbld h936iJ) and lZwickvl 
Traditionally the favoured mechanisms have been eyeball 
classification against a set of visual criter ia 1 Hubble 19 26), re- 
sulting in the [Hubble] tuning fork jjeansi[l92 9. Hubbl elll936i^ 
Sandagelfl96ll) . or similar schemes (e.g.. Ide Vaucouleurg |l956 



eil) , 

van den Berghlll97dj. and via lig ht-profile fi tting ( de Va ucouleurd 
194^. 1959: 'S ersidll96l 1968; 'Freeman T970). Both systems 
remain in common usage and are now being routinely applied 
to large datasets in an automated manner. For example. Artifi- 
cial Neu ral Networks are now used to im plement eyeball c lassi- 
fication <Naim et alJll995l: IOdewahrJll9 95: Ball etajl Eooi) and 

i iublic co des such as GIM2D (Simardetal. 2002) and GALFIT 
Peng et al. 2002) are now available for automated bulge-disc de- 
composition or single Sersic profile fitting. Much of the drive to 
pursue these methods now comes from our ability to apply coarse 
structural measurements at any redshi ft via high-resolutio n space- 
based observatories (see for ex ample D river et alJll995S 1995b; 
IlIuv et al.. 1998: Ravindranath et aL2004i) ! 

Significant effort is also being invested in exploring al- 
ternatives such as: the joint luminosity-surface brightness(size) 
plane (LSP; see JJriver et al. 2005 and references therein), the 
colou r -luminosity plane (CLP; Strateva et al. i200lt iBaldrv et alj 
|2004 iBalogh et all 120043: Ipaljer et al. 2005) — both of which 
have traditional roots^ — spect ral classificati on iMadgwick et alJ 
200?), line strengths I' Linetal.i il996i: iKauffmann et alJ l2003l 
Balogh et al. 2004b), Concentration/Asymmetry/Smoothness 

J iarameters jC onselice et al. l2004h. the Gini coefficient 
Abraham, van den Bergh & NaiJ l2003l iLotz. Primack & Madai] 
|2004 , Pri nciple Component Anal vsis iElUs et alj|2005h . Shapelet 
analysis I'Kellv & McKaJ l2005l) . and Fourier decomposition 
(Odewahn et al. 2002) for example. At present no single approach 
stands above the rest and all are open to various criticisms. For 
example eyeball morphologies are subjective, the LSP and CLP are 
too coarse, Concentration/Asymmetry/Smoothness is susceptible 
to short term transitory effects (minor mergers, interactions, 
episodic star-burst etc), the Gini coefficient and PCA analysis are 
unlikely to lead to a straightforward connection to the underlying 
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physics, not all galaxies are readily profiled by Sersic or Sersic 
plus exponential models, and Fourier decomposition requires 
high signal-to-noise and is hence not applicable beyond the local 
universe. In addition no method has been shown to have a deep 
phy sics basis, a lth ough galaxy profiling p robably comes closest 
(see lKinJ 19631 and lFall & EfstathioJl98cl for example). 

Moving beyond classification issues, the standard method 
for representing the collective galaxy population has typically 
been via the luminosity distribution (e.g.. .Hubble. 1 936h: Binggeli, 
Sandage & Tammann 1988) and its adop ted analytical rep resen- 
tation, the Schechter luminosity function ^Schechtei"l976'). Con- 
temporary measurements of the global B-band galaxy lumi- 
nosity distribution started with the CfA slice 1 Davis & Huch^ 
fl982i) followed by the analysis of lEfstathio~et aL L198 8) and 
have culminated in recent results from the ESO Slice project 
iZucca et al. 1997), the Two-degree Field Galaxy Redshift Survey 
(Norbera et al. 2002), the Sloan Digital Sky Sur vev l^Blanton et alj 
2003b), and the Millennium Galaxy Catalogue ( iDriver et al 120051) 
and references therein. Generally these surveys now concur on the 
space-density of luminous galaxies (Mbmgc ~ 51og/i < —16 
mag), with the space-density of dwarf systems {M s^^c ^ ^ l^S ^ > 
— 16 mag) essentia lly unconstrained (Driver 2004, although see 
iBlanton et all20o3 for a considered attempt based on SDSS data). 
Data from each of these surveys has also been used to provide 
galaxy luminosit y distributions subdiv ided by various criteria, e.g ., 
spectral type ( Madgwick et al. 2002), colour iZuccaetal. 199"^, 
galax y light-profile sha pe (Blanton et al. 20-03b.) and morphology 
iNak amura et al. ll2003^ . These and earlier attempts to subdivide 
the galaxy population have been comprehensively summarised by 
^e Lapoarent (2003), who highlights the complexity and confusion 
that can arise through the comparison of galaxy LFs subdivided by 
differing criteria (i.e., although morphology, spectral type, colour 
and profile shape have long been known to be correlated, the terms 
are not synonymous). This once again leads back to the question 
as to what is the fundamental way to divide the galaxy popula- 
tion if at all? This is not an easy question to answer but empiri- 
cally one can address it by highlighting the division that leads to 
maximum variance in the recovered Schechter function parameters 
as well as by highlighting modality in joint distributions. Recent 
studies have identified bimodalities in the colour-magnitude plane 
( Strateva et al. 2001, Ball et al. 2005), the stellar age-stellar mass 
plane tKau&iani^^alj|2002|), the luminosity-surface brightness 
plane jPr^e^^HoOg) and the luminosity-size plane iShen et all 
I2OOI 'Blanton et al.'2003aV While the nature of these bimodalities 
is no doubt of common origin it remains unclear as to which plane 
is more fundamental and what these bimodalities might be telling 
us. Our argument, which we develop throughout this paper, is that 
they all reflect the two component nature of galaxies (as opposed to 
two distinct classes of galaxies) and that the fundamental division 
is between the disc and bulge components, each of which may have 
a distinct formation mechanism acting over two distinct eras. 

U sing the Millennium Galaxy Catalogue (MGC: lLiskeetalJ 
l2003t) we explore the coarse global properties of luminous galax- 
ies {Mbuqc — 51og/i < —16 mag). In Section 2 we summa- 
rize the MGC and describe the morphological classification pro- 
cess which supplements the colour, spectral and struc tural param- 
eters which are described in detail elsewhere (see iDriver et alJ 
[2005, Allen et al. 2005). In Section 3 we explore the distribution 
of galaxy morphological type, MGC continuum type, rest-(it — r)- 
colour, global Sersic profile shape and derive the luminosity distri- 
butions divided along natural boundaries. In Section 4 we continue 
the exploration by examining selected bivariate distributions. Im- 
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plications for galaxy formation are discussed in Section 5. We use 
Q.M = 0.3, r^A = 0.7, Ho = lOO/i km s"^ Mpc~^ throughout. 



2 THE MILLENNIUM GALAXY CATALOGUE 

The Millennium Galaxy Catalogue fMGCi lLiske et aljEoOsI) is a 
deep (/^LiMiT = 26.0 -Bmgc mag arcsec"^) survey of a 37.5 deg^ 
region of sky, ~ 0.6 deg wide and extending from ~ 10 to 
~ 14'' 50™ along the J2000.0 equator. The imaging survey was 
conducted with the Wide Field Camera installed on the 2.5m Isaac 
Newton Telescope in La Palma and the survey's design, execu- 
tion, redu ction, object detec tion, and preliminary analysis are de- 
scribed in lLiskeetalj|2003h . The MGC lies within the Two Degree 
Field Galaxy Redshift Survey (2dFGRS) Northern Galactic Cap 
region and the Sloan Digital Sky (SDSS) Early Data Release Re- 
gion. A detailed comparison of the MGC w ith the much larger bu t 
shallower SDS S (Data Release On e, DRl. lAbazaiian et alfcoosi) 
is described in ICross et all i2004) . The spectroscopic extension 
(MGCz), which built upo n the redshift s provi ded by the 2dFGRS 
and SDSS, is described in lDriver et alJ i2005l) . The MGC-BRIGHT 
imaging catalogue contains 10 095 galaxies (13.0 < -Bmgc < 20.0 
mag), for which redshifts have now been obtained for 9 696 result- 
ing in a global completeness of 96 per cent. In this paper we oc- 
casionally restrict ourselves to galaxies brighter than Bmgc < 19.0 
mag. In these cases the sample size is reduced to 3 492 galaxies for 
which redshifts have been obtained for 3 487 indicating a global 
completeness of 99.9 per cent. Selection biases wit hin the cata- 
logue and their treatment are extensively discussed in lDriver et alj 
t2005,) . The MGC photometric system (-Bm gc) is calibra t ed to V ega 
and common filter conversions are shown in iLiskeetalJliool and 
Appendix A o f ICross etalJj2004 t). The MGC catalogue used in this 
paper is available online^ and includes 120 parameters per galaxy 
derived from our own software, GIM2D and through matching to 
the 2dFGRS and SDSS-DRl datasets. We now describe the param- 
eters used in this paper and their derivations starting with eyeball 
morphology. The resulting observed distributions (i.e., uncorrected 
for volume bias) of each parameter (morphology, 2dFGRS rj, MGC 
continuum classification, global and core colours, mass-to-light ra- 
tio, Sersic index, central surface brightness and half-light radius), 
are shown on Figs.Qand|2|for the -Bmgc < 20, A'/bmgc ~ 5 log h < 
— 16 mag and -Bmgc < 19, Mbmgc ~ 5 log /i < — 16 mag samples 
respectively. 



2.1 The MGC morphologies 

We have derived the morphological classifications for the 3 314 
galaxies that lie within the MGC si ze and surface brightness se- 
lection boundaries (see lDriver et alJt200 5) with -Bmgc < 19 mag 
through eyeball classification (SPD). This is a subjective process 
and was achieved by visual inspection of batches of 20 grey scale 
images displayed over three surface brightness ranges: 19.0 - 26 
mag arcsec~^, 20.5 - 26 mag arcsec^^, 22.0 - 26 mag arcsec^^. 
The three visual classes are defined as follows: 

• E/SO: Smooth, highly concentrated symmetrical systems with no 
obvious spiral-arm disc component. 



^ http://www.eso.org/~jliske/mgc/ 



• Sabc: Clear two-component bulge-disc system, generally smooth 
and symmetrical. 

• Sd/Irr/Pec: Disc only (i.e., no obvious single core), asymmetric, 
highly disturbed or multiple-core systems. 

These classifications are consistent with our deeper Hubble Space 
Telescope studies of more distant galaxies ( Driver et al. 1995a|, 
1995b, 1998, 2003: .Driver 1999; Odewahn et al, 1996 : Cohen et aU 
2003), which are based on comparable physical resolution and 
signal-to-noise data. The entire classification process was under- 
taken by SPD over a period of some weeks and repeated several 
times until no reclassifications were necessary. Fig. |3] sho ws rep- 
resen tative examples of each class including stars (see .Liske et alJ 
l2003l for details of the star-galaxy separation). Note that images of 
both red and blue 'ellipticals' are shown and these are discussed 
in Section 3.2. Table Q contains the morphological number counts 
which can be compared to our earlier and deeper HST studies listed 
above. 

2.2 The 2dFGRS rj parameter 

The 2dFGRS 77 parameter is defined in lMadswick et al] i2002h and 
is a linear combination of PCI and PC2 from a Principle Com- 
ponent Analysis of the spectra. PCI contains information on the 
emission and absorption line strengths and the continuum (with 
roughly equal weights) and PC2 contains information on just the 
line strengths. The linear combination maximises line features 
and hence provides a good indicator of current and previous star- 
formation. The reader is referred to fig. 7 of Madgwick et alJ^2002^ 
for canonical examples of rj type 1, 2, 3, & 4 spectra. 

2.3 MGC continuum classification 

We divide the sample into three MGC continuum types using our 
spectral template fitting code (see Driver et al. 2005) and restricting 
it to fit either a 15 Gyr elliptical (El), a 7.4 Gyr early-type spiral (Sa) 
or a 2.2 Gyr late-type spiral (Sc). All three continuum spectra are 
taken from Pogaianti 1 1997). Fig.|4|shows both the location of the 
broad band filters and our three canonical continuum shape spectra 
shifted to 2; = 0.1. The op timal fitting s pectru m for each galaxy 
is obtained as described in lDriver et alj lioosi) . In brief, the flux 
for the uBgriz filter set is derived for each template at the known 
redshift of the object and ^^-minimisation used to determine the 
best fitting template (where the normalisation is marginalised over). 
Where no satisfactory fit is found the filters are dropped in the or- 
der of -Bmgc (non-SDSS), then u (red leak) and then z (fringing). 
The optimal spectra provide an indication of the continuum shape. 
In essence the MGC continuum typing is simply a sophisticated 
colour cut whioch uses all available colours hence correlation be- 
tween continuum type and colours is to be expected. 

2.4 SDSS global and core colours 

Overlap with the SDSS-DRl enables complim entary uqriz pho - 
tometry for 99% of the MGC galaxy sample (see'Cross et al.''2004|). 
To convert these colours to rest-wavelengths we use the appropri- 
at e spectral template and a simple evolutionary recipe as described 
in'Driver et all j2005l). The SDSS -DRl colour bimodality was first 
reported bvTstrateya et alJ i200lh and studied in more detail by 
ISaldrv et alj J2004 . Its environmental dependency was explored 
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Figure 1. The various criteria by which the MGC can be subdivided (as indicated, see also Sections 2.1 - 2.8 ). The shading represents the location of MGC 
galaxies classified based on their continuum shape as either El (solid), Sa (hashed) or Sc (diagonal lines). 



further in lBaloghetalJj2004j) . These studies typically adopted the 
global {u — r) rest colour. Here we show both the global and the 
core (u — r) rest colours (hereafter {u — r)g and {u — r)c re- 
spectively). The former is derived from Petrosian magnitudes and 
the latter from PSF magnitudes. One may argue whether the core- 
colour is better defined by the PSF magnitude, by a minimal fixed 
aperture colour, or even fibre magniutdes which are degraded to 
a consistent poor seeing value. Clearly, PSF magnitudes are more 
appropriate for ultra-compact sources, and a fixed aperture slightly 
larger than the seeing more appropriate for extended sources. How- 
ever as the cores are often barely resolved we essentially opt for the 



smallest credible measurement which is the PSF magnitude. While 
this is somewhat profile dependent an equivalently small fixed aper- 
ture would be biased by the wavelength dependency of the PSF 
size. Both of the colour distributions are seen to be strongly bi- 
modal (see Fig.Ql and^, respectively) with the peaks marginally 
better separated in the {u — r)c distribution. We separate the galaxy 
population into two classes with cuts at {u — r)g = 2.1 and 
(u — r)c = 2.35 respectively. 
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Figure 2. As for Fig.C]except for galaxies with Bmgc < 19 mag only throughout. The shading now highlights how the eyebal morphological classifications 
into either E/SOs (solid), Sabcs (hashed) or Sd/IiTS (diagonal lines). 



2.5 Stellar mass-to-light ratios 

iBell & de Jond i200lh provide a standard prescription to deter- 
mine stellar masses and stellar mass-to-light ratios from broad band 



colours. From their Table 1 we have: 



XogiM' /Lb) = -1.224 + 1.251(B - R). 



From lCross et alj i2004) we have the following conversions from 
SDSS-DRl photometry of: 

(B -R) = 1.54(g -r) + 0.35, 



and can recast the above as: 



(A'/7Lb) = 10[^-^3(''-'-'-o-^^1, 



where M* / Lb is the -Bmgc stellar mass-to-light ratio in solar units. 
We elected to use the SDSS-DRl colours for consistency rather 
than (-Bmgc — r) because the differing deblending algorithms occa- 
sionally lead to spurious (Bmgc — r) colours. Not surprisingly, as 
the stellar mass-to-light ratios are derived directly from the global- 
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Figure 3. A random sample of our postage stamp images for representative examples of (left to riglit) stars, red ellipticals, blue ellipticals (dEs?), spirals, 
and late-type irregulars at various magnitude intervals (as indicated). Faintwards of Bmgc = 19 mag morphological classification becomes intractable and 
higher-resolution data in better seeing is required. 
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Table 1. Morphological number-counts from the MGC. 



-Bmqc (mag) log N[A11] log N[E/SO] log N[Sabc] log N[Sd/lrr] 



13.75 
14.25 
14.75 
15.25 
15.75 
16.25 
16.75 
17.25 
17.75 
18.25 
18.75 
19.25 
19.75 



-1.012 
-0.790 
-0.489 
-0.343 
0.123 
0.396 
0.677 
0.894 
1.228 
1.474 
1.686 
1.901 
2.122 



-1.188 
-1.489 
-1.012 
-0.790 
-0.535 
-0.127 
0.226 
0.385 
0.717 
0.978 
1.213 



-1.489 
-1.012 
-0.644 
-0.586 
-0.074 
0.200 
0.373 
0.640 
0.911 
1.143 
1.357 



-1.489 

-1.489 
-0.711 
-0.790 
-0.147 
0.015 
0.547 
0.802 
0.974 
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Figure 4. The three adopted continuum spectra used in our fitting process. 
The spectra are shifted to a redshift of 0.1 and are arbitrarily scaled. We 
also show the transmission curves of the filters for which photometry is 
available. 



the distribution of the best-fitting Sersic parameters for the MGC. 
The distribution in log(n) is bimodal, presumably reflecting bulge 
dominated (i.e., pressure supported) or disc dominated (i.e., rota- 
tionally supported) systems. Note that it is more appropriate to plot 
log(n) rather than n, as n appears in the exponent of the intensity 
equation, i.e., a small change in n makes a large difference for low 
n values but minimal difference for high-n values. We adopt a di- 
vision at the saddle point of n = 2.0. From repeat observations 
of 702 objects we have quantified the random error on the Sersic 
index to be A log(n) = ±0.041. 

2.7 The MGC central surface brightnesses 

The (extrapolated) central surface brightness (^o) is calculated di- 
rectly from the GIM2D Sersic only fits and hence strongly de- 
pendent (although not exclusively) on the Sersic index. It is in- 
cluded here for completeness as many earlier studies have subdi- 
vided the population by a surface brightness cut. The distribution 
is bimodal and we divide the galaxy population at = 19 i?Moc 
mag arcsec"^. Note that in our previous studies ( ICross et alj200l|: 
fCross 8l Driver 2002; Driver et al. 2005) we explored the effective 
surface brightness absolute magnitude plane. The distinction is im- 
portant, whereas bulges and discs typically have similar effective 
surface brightnesses their central surface brightnesses are signifi- 
cantly different. For overcoming selection bias the former is more 
relevant, for exploring variations the latter is more appropriate. 



(cj — r) colour, we again see (Fig.Q^ a bimodality with an obvious 
division atlog(M*/is) = 0.45. 

2.6 The MGC Sersic indices 

The MGC Se rsic profile indices (n), are derived using GIM2D 
fSimard et al.' 2002) and the process is described in detail in 
kllen et al- C200SV Briefly, GIM2D performs a x^-minimisation 
between the image and a 6-parameter model image. A segmenta- 
tion mask derived from SExtractor is used to indicate which pixels 
are associated with the object, and an appropriate PSF image is 
convolved wit h the model i mage. The software use s the standard 
Sersic profile ISersidll96ilGraham & Driverll2005ft returning the 
total flux, Sersic index (n), half-light radius (7?e), ellipticity, po- 
sition angle and positional and background offsets. Fig.^ shows 



2.8 The MGC half-light radii 

The half-light radius {Re), is calculated directly by the GIM2D 
code and is therefore PSF-corrected. The distribution is relatively 
broad and peaked at around 10 kpc. Hence, while bimodalities exist 
in many observable properties, the size distribution is smooth (and 
hence the effective surface brightness as well). There is no obvious 
location at which to divide the population. 

2.9 The observed (uncorrected) distributions 

Fig. Q shows the observed distributions, i.e., uncorrected for vol- 
ume bias, for galaxies with Bmgc < 20 mag and Mbmgc ^ 
5 log ft < —16 mag, except for morphological type and r\ which 
are only known to Bmgc = 19 mag. Note that the absolute magni- 
tude limit essentially confines this study to the giant galaxy regime 
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(with some contribution from thie brightest dwarf ellipticals). In 
Fig. |2| the overall distributions also shown subdivided according 
to MGC continuum type (solid. El; hashed, Sa; or diagonal lines, 
Sc). This highlights the correlated nature of the bimodalities. How- 
ever, it is also clear that the subdivisions of the galaxy population 
based on colour, morphological type, etc., will each extract distinct 
subsets of the galaxy population. Hence, great care must be taken 
when comparing high and low redshift samples subdivided by dif- 
ferent methods. Fig.|2|shows the same distributions but restricted to 
the Bmqc < 19 mag population. The morphological groups (E/SO, 
Sabc, Sd/Irr) are shown with distinctive shading throughout (solid, 
hashed, diagonal lines, respectively). From Fig.|2|it is clear that the 
E/SO (bulge-dominated) and Sd/Irr (disc-only) populations occupy 
the extremes of the bimodal distributions with the Sabc distribution 
as a bridging population. 



3 THE LUMINOSITY DISTRIBUTIONS AND 
SCHECHTER FUNCTIONS 

To obtain luminosity distributions for sub-populations we first re- 
derive the global bivariate brightness distributions for the i?MGc < 
19 ma g and Bmgc < 20 mag sa mples following [ Driver et alJ 
i2005h . The method implemented in lDriver et alj i20()5l) is a vari- 
ant of the standard s tep-wise maximum likelihood (SWML; see 
lEfstathiou et all 1988 ) operating in the plane of luminosity and sur- 
face brightness (as opposed to luminosity alone). The main ad- 
vant age is the ab ility to track the five key selection boundaries 
(see lDrivej|l99 ^ which govern whether a galaxy is reliably de- 
tected or not: maximum/minimum luminosity, maximum/minimum 
size, and maximum mean effective surface brightness. The anal- 
ysis'' assumes (^m = 0.3, SIa = 0.7) with Ho = lOOh km 
s~^ Mpc~^ and K-corre ctions derived individually for each galaxy 
(see iDriver et a l."2Q05' for specific details). Fig. |^ shows a com- 
parison of the luminosity distributions, functions, error ellipses and 
LSP distributions for these two samples which are fully consistent 
with each other. 

To derive the subdivided luminosity distributions, tp{M, x) 
we now simply calculate: 

fi — oo 

^{M,x)= J2 [MM,^i).<j,{M,^^)]Sfi, (1) 

where fx {M, represents the apparent frequency of some at- 
tribute, X, (e.g., fraction of E/SO type) within that absolute mag- 
nitude and surface brightness bin, SMSfi. 4>{M,y) represents the 
space density of all galaxies in that same interval. The integral 
over surface brightness then converts the bivariate distribution to 
the monovariate luminosity function; note that although the inte- 
gral ranges across all values of surface brightness the distribution 
at any M is actually bounded by the selection limits. This is essen- 
tially identical to re-deriving the LSP (also known as the bivariate 

^ In detail a few minor changes have occurred in the analysis since 
iDriver et alJ l2005t) . Firstly the selection boundaries are implemented he- 
fore the seeing correction; secondly the incompleteness matrix is also 
derived prior to seeing correction; and thirdly the normalisation is now 
made to galaxies in the range: —20.0 < Mb^hqq < —19-5 mag and 
20.75 < (Atgff) < 23.25 mag arcsec'^and 0.047 < z < 0.162. Minor 
improvements to the code have also resulted in a more accurate evaluation 
of LSP cells which are intersected by the selection limits. These updates do 
not change the results o flDriver et alj <2005l) significantly. 



brightness distribution or BBD, see lBovce&PhilliDPslll995h for 
each sub-population but computationally simpler and more stable. 
Figs. |6||9| show the recovered luminosity distributions for the var- 
ious subdivisions. In each case the luminosity distributions were 
fitted by a Schechter function to Mbkgc " 5 log h = —16 mag. In 
each figure the recovered data are shown (left) along with the best 
Schechter function fits and the l-,2- and 3- a error contours (right- 
side). The fitting parameters are also tabulated in Table |2| along 
with the 1-(T errors from the fit, the integrated luminosity den- 
sity (jb ), the median (it — r)g and the total stellar mass for each 
subdivision. The luminosity densities and stellar densities (M*) 
were derived as follows: 

jb, =</>*io-°-''(^^^-*^®'r(a + i), 

and, 

j^/* ^ ]^Q[1.93(9-r)-O.79)-j^g-O.4{MB-M0)^ 

i.e., analytically for the luminosity density and an empirical sum for 
the stellar mass density because of the strong colour dependency. 

3.1 Luminosity distributions/functions by morpliological 
type 

In Fig. |S| (upper) we show the eyeball morphological luminosity 
distributions/functions. Although the three divisions show segrega- 
tion in their Schechter function parameters, the distinction between 
the E/SO and Sabc categories is not particularly strong. Certainly 
the E/SOs do not appear to follo w the Gaussian distributi on an- 
ticipated bv lBinggeli et alJ Il98^ . r eported in Jerien & Ta mmannI 
1 1997), and discussed extensively in lde Lapparenl^ 1.2003.,) . In Ellis 
et al. (2005, their fig. 3) we found that our eye was often fooled by 
smooth, blue, low-luminosity spheroidal systems. This highlighted 
both the subjective nature of eyeball classification and its suscepti- 
bility to variations in data quality. Whether the blue spheroid pop- 
ulation constitutes the tail-end of the 'downsizing' phenomenon or 
interlopers is at this stage unclear but under investigation. If one 
uses additional colour information to separate out the blue and red 
spheroids (i.e., following Ellis et al.) and re-derives the morpho- 
logical luminosity distributions we now get the distributions shown 
in Fig.|6|(lower). The luminosity distribution of the blue spheroid 
(BS) population alone (not shown but see Table |2} follows almost 
exactly the same distribution shape as the Sd/Irr population and 
most likely has a common origin. This is also consistent with the 
idea that the E/SO population has indeed been contaminated by 
a population of smooth blue dwarf systems (dwarf spheroidals). 
In so far as a comparison can be made, our valu es could be con- 
trasted to those listed in Table 4 of Ide Lapparenlf 1 2003) where we 
see measurements for the E/SO class from t he Stromlo APM sur - 
vey (A/smgc = -19.33 mag, a = +0.2; iLovedav et alJll992l) . 
the Southern Sky Redshift Survey 2 (M%^^^ = -19.64 mag, 
a — —1.00; Marzke et al. 1998) and the Nearby Optical Galaxy 
Sample iMg^^^ = -19.80 mag, a = -0.97; iMarinoni et alJ 
119991) . These three surveys concur well with our uncontaminated 
and contaminated E/SO samples suggesting the problem of BS con- 
tamination may be endemic in all eyeball classified samples (and 
hence in ANN classification as well) - see E llis et al. 1 2005) for 
the formal detection of these systems and also lCross et alj i2004h 
for discussion of BS contamination of field ellipticals at z ~ 0.75 
(r.f., fig. 2a from Bell et al. 2005). While arguably a combined eye 
and colour cut may provide a more robust division, the process of 
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Figure 5. (main panel) the global B-band luminosity function derived via our modified step-wise maximum likelihood method for the Bmgc < 19 (solid 
squares and solid line) and -Bmgc < 20 (open circles and dotted line) mag samples, {upper right) the 1,2 and 3 a error-ellipses of the Schechter function fits 
and {lower right) the luminosity surface brightness plane and the logarithmic contours for the Bmgc < 19 (dotted line) and Bmgc < 20 mag (solid line) 
samples. The two dashed fines denote the limits of reliability for the two samples. 



morphological classification (whether by eye or ANN) is subjective 
and hence susceptible to confusion, controversy and ambiguity. 

3.2 Luminosity distributions/functions by r/ and MGC 
continuum type 

Fig. shows the luminosity distributions/functions subdivided by 
T) (upper) or MGC continuum type (lower). The luminosity distri- 
butions/functions subdivided b y r? type concur with the earlier re- 
sults by the 2dFGRS team (i.e.. lFolkes et alJl999l:lMadgwick et alj 



|20o3), with the ryl-type showing a distinct luminosity distribution 
and less distinction between the remaining 77 types. Unfortunately 
the fraction of galaxies without a known rj type is significant (even 
to B < 19 mag) due to the requirement for moderately good signal- 
to-noise spectra. 

Our continuum typing, based on spectral template fitting to 
the available uBgriz broad-band photometry, appears significantly 
more promising with the three continuum shapes (see Fig.|4} rep- 
resented by very distinct luminosity distributions/functions. From 
Fig.Qwe see that these distinctions (shown as shaded, hashed and 
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Figure 6. (upper) The derived data and fitted Schechter functions subdivided by morphological type and (lower) with the blue spheroids (BS) moved into the 
Sd/Irr categoiy. 



lines histograms) are well correlated with other observables. While 
correlation between continuum type and colour (and hence stellar 
mass-to-light) is to be expected, the correlation with structural mea- 
surements is significant as no structural information was used in 
the determination of the MGC continuum type. The El continuum 
type extracts the red, high mass-to-light ratio, concentrated, and 
predominantly luminous systems. The Sa class exhibits a broader 
range bridging the El and Sc populations. The Sc class constitutes 
the bluest and least concentrated systems. The implication is that 
the El systems are luminous (massive) and old (or at least their stars 



are), and the Sc systems appear to be predominantly low luminosity 
(low mass), and young (late-formers). 



3.3 Luminosity distributions/functions by colour and 
mass-to-light ratio 

Fig. m shows the luminosity distributions/functions subdivided by 
{u — r)g (upper) and (it — r)c colours (lower). It is worth noting 
that our calculation of the rest colour is independent of the standard 
SDSS method, based on our own spectral temp late fitting routin e 
which includes all 27 spectral templates given bv lPoggiantil il997h . 
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Figure 7. {upper) The derived data and fitted Schechter functions subdivided by 2dFGRS-T7 parameter and {lower) subdivided by MGC continuum type. 



see lPriver et alJ i2005 ^ for details. The luminosity distributions of 
the red and blue populations are significantly distinct. Worth noting 
is that the (u — r)^ colour - derived from SDSS-DRl PSF magni- 
tudes - appears to segregate galaxies better than the global colour - 
derived from SDSS-DRl Petrosian magnitudes. The improvement 
in segregation (see Fig.|8|right side panels) is significant and could 
be interpreted as better separating the galaxy population into those 
with dominating old bulges and those without. The key question is 
why do red bulges predominantly occur in luminous galaxies? Pre- 
sumably either the low-luminosity bulges are unresolved and hence 
hidden within the PSF, or that old inert (and therefore red) low- 



luminosity bulges are rare*. This is consistent with some studies 
which argue that the bulges of late-type systems are ps eudo-bulges 
built up through inner disc and bar instabilities (e.g.. lErwin et alj 
■2003.) rather than genuine bulges. As genuine bulges should be cus- 
pier one might expect the core colour to be better at identifying such 
systems. As mentioned in Section 2.4 our core colour is not ideal 
but simply relies on the colour within the smallest available aper- 
ture. An obvious extension would be to derive the core colour from 



* Although dwarf ellipticals are considered the natural extension of the 
disc-less spheroids, Graham & Guzr nan 2003 , there is no obvious low- 
luminosity bulge-disc counterpart. 
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Figure 8. The derived data and fitted Schechter functions subdivided by (w — r)g colour {upper) and (u — r)c colour {lower) 



multicolour 2D bulge-disc decomposition models to see whether 
the red peak becomes even sharper. 

The SPSS iu — r )a colour bimodality was studied in detail 
by iBaldrv et alj i2004) who derived the r-band luminosity func- 
tions and identified faint-end parameters of a = —1.35 ± 0.05 
and a — —0.83 ± 0.002 for the blue and red populations re- 
spectively. This compares with our values of a = —1.35 ± 0.02 
and a — —0.54 ± 0.05 for our B-band selected samples. As 
one expects a _B-band sample to preferentially miss red systems 
when compared to an r-band sample, the higher a for the red 
systems is unlikely to be significant. As the mas s-to-light ratios 
are derived from the global {g — r) colour ("via iBell & de Jon3 



I2OOII) the narrowness of the mass-to-light ratio distribution is ex- 
pected. However, the MGC continuum classification appears re- 
markably robust at identifying systems occupying the high mass- 
to-light ratio peak (Fig. Q^. Perhaps even more interesting is the 
extreme narrowness of the high-M/L peak, essentially a represen- 
tation of the well known colour-magnitude relation of cluster el- 
lipticals iSandage & Visvanathan 1978). This narrowness presum- 
ably reflects the universal endpoint of galaxy evolution, which in 
turn places a joint constrai nt on the age, metallicity and even the 
universality of the IMF (see lBruzual & Charloi2003h . 
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3.4 Luminosity distributions/functions by structure 

Fig. |9| shows the luminosity distributions/functions subdivided by 
Sersic index, n (upper), and central surface brightness (lower). As 
the central surface brightness is derivable from n, total magnitude 
and effective radius, it is not independent and hence the similarity 
in the two distributions (upper and lower panels) is to be expected. 
Taking the log(n) division we see that while the high and low-n 
populations follow distinct luminosity distributions, the distinction 
is not as great as for the MGC continuum and colour divided distri- 
butions. This is, at least in part, due to the fact that some dEs with 
— 16 > Mbucc ~ 5 log /i > — 18 mag will have n < 2 and others 
will have n > 2. Folding in colour information to separate the red 
dEs from the blue discs will help and is considered in the following 
section. 



3.5 Luminosity distribution summary 

The main conclusions from this section are: 

(1) Morphological classification based on the eye, including Arti- 
ficial Neural Networks which use a training set, are fraught with 
difficulty because of the subjective nature of the classification pro- 
cess. 

(2) In so far as morphological luminosity functions can be obtained, 
particular care must be taken to distinguish genuine ellipticals from 
smooth, blue, low-luminosity systems and eyeball classifications 
should be made using colour images if possible. 

(3) Both global and core (it — r)g,c colours segregate the population 
well into blue and predominantly luminous red systems. 

(4) The core colour appears to better segregate the two populations 
(in M — a-space) than the global colour, probably due to the blend- 
ing of colours from red bulges embedded in blue discs. 

(5) The narrowness of the red, core-colour peak may imply an early 
formation age from a universal initial mass function for these sys- 
tems (as commonly assumed, e.g.. Bell & de Jons 2001). 

(6) Further gains may be made by moving beyond a single colour 
cut to a continuum-type classification system based on multiple 
colours. This highlights three distinct populations. A luminous 
red population a low luminosity blue population and an interme- 
diate population, these are closely associated with E/SO (bulge- 
dominated), Sd/Irr (disc-only) and Sabc (bulge and disc) systems. 

(7) Galaxies can also be effectively separated into two populations 
using structural division in log(n) space. The populations consti- 
tute a concentrated and diffuse population. 

(8) Colour appears to segregate the galaxy population more effec- 
tively than any other single measurement (i.e., Sersic index or sur- 
face brightness). 

In the next section we continue the analysis by exploring 
selected bivariate distributions. 



4 BIVARIATE DISTRIBUTIONS 

We derive volume corrected bivariate distributions/functions in a 
similar manner to the luminosity distributions by starting from the 
overall bivariate brightness distributions for the entire Bmgc < 19 
and Bmgc < 20 mag population and projecting attributes x and y 
on to a new plane, i.e.. 



where the values are defined as for Eqn. (1). Fig. ll0tIT4l show the bi- 
variate distributions of: {u — r)g, {u — r)c, log(n), and log(_Re) 
versus M b^gc ■ ^'"-'^ figure we show (upper left) the raw ob- 
served distribution as both data points and contours, and (lower 
left) the volume corrected distributions using Eqn. (2). The right- 
side panels show the histograms along the absolute magnitude axis. 
As a reminder, only galaxies with Mbmgc ~ 51og/i < — 16 mag 
are included in the analysis. 

4.1 Bivariate distributions by colour 

Figs.^|and II ll show the bivariate distributions of rest {u — r)g and 
(u — r)c colours versus Mb^gc- both figures we see a dis tinct 
red and blue sequence as characterised bv'Bald ry et al.l J2004) us- 
ing SDSS data. In both figures the red distribution appears entirely 
bounded albeit overlapping with the luminous end of the more ex- 
tensive blue population. Both the red and blue sequences show 
colour gradients reminiscent of the luminosity-metallicity relation 
seen in nearby rich clusters. Note that bimodality is only obvious 
in the observed histograms and not in the volume corrected his- 
tograms as the numbers are dominated by the lower luminosity sys- 
tems. In terms of stellar mass (not shown) the bimodality remains 
but is skewed towards the red peak. The division between the red 
and blue populations is more distinct (and the red peak histogram 
narrower) when (u — r)c colour is used. This implies that the core 
colour is the more fundamental (less contaminated) measure. Fur- 
ther data probing to lower luminosities is required to unequivocally 
conclude that the red sequence is bounded at low luminosity (i.e., 
few galaxies with {u — r)c > 2 and AIbugc " 51og/i < —17 
mag). An indication of a low luminosity red spike (dEs ?) is seen 
but with extremely low significance because of limited statistics 
and the large volume amplification at the faint-end. 

4.2 Bivariate distributions by structure 

Figs. ll2lll3l and fT4l show the corresponding structural bivariate dis- 
tributions for log(n), Ho and i?e respectively. The first two show 
bimodality albeit less obviously than for the global and core-colour 
distributions. Both Figs. 1121 and 1131 suggest a concentrated, high 
surface brightness, low-luminosity population close to the limit of 
analysis (in the volume corrected distributions). Whether this rep- 
resents a third distinct population or a manifestation of the limiting 
statics close to the selection boundaries is difficult to establish and 
requires additional data to clarify. One clear possibility is that it 
may represent the onset of nucleated dwarf systems which would 
exhibit both higher than typical Sersic indices and correspondingly 
higher central surface brightnesses. Unlike the colour bivariate dis- 
tributions, the distinct concentrated giant (E) distribution blends 
into the diffuse (low-n, disc) population towards lower luminosi- 
ties. The size distribution shows no bimodali ty but a bro aden- 
ing towards lower luminosity as found by .Shen et alj J2003j|) and 
lDriveretalJj2005h . 

4.3 The joint distribution of core-colour and Sersic index 

Having established that bimodality exists in both the colour and 
log(n) distributions, and that the {u — r)c shows the bimodality 
most strongly, we now explore the joint (it — r)c-log(n) plane. 
Fig. 115! shows the observed distribution (upper) in contour form 
(left) and as a 3D projection (right). The lower panels show the 
equivalent for the volume corrected distribution. Fig. ll5l shows two 
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Figure 9. The derived data and fitted Schechter functions subdivided by log(Sersic index) (upper) and fitted central surface brightness (lower). 



distinct peaks with a bridge between them pointing towards either 
two distinct populations or two dominant characteristics. The peaks 
are typified by a high Sersic index red peak and a low Sersic index 
blue peak. Note that the additional third spike (in the lower right 
panel) may constitute a compact blue population reminiscent of 
the blue spheroids (see Section 3.1). The key question at this stage 
is whether we are seeing two distinct galaxy populations (red and 
blue) or whether the peaks correspond to two characteristics (i.e., 
bulges and discs). Bulge-disc decomposition is clearly essential 
to establish this unambiguously. However early indication comes 
from the morphological classifications where, by definition (see 
Section 2.1), E/SOs represent bulge-dominated systems, Sabcs joint 



bulge-disc systems and Sd/lrrs disc-only systems. Fig. 1161 (left) 
shows the raw observed distributions (i.e., equivalent to Fig. llSi rup- 
per right]) for the E/SO (upper), Sabc (middle) and Sd/Irr (lower) 
populations. We see quite clearly that the bulge-dominated (E/SO) 
and disc only (Sc) systems lie almost exclusively in one peak or 
the other. The bulge-l-disc systems (Sabcs) straddle both peaks and 
the divide. This is extremely strong evidence that points towards 
the fundamental nature of bulges and discs and the need to inter- 
pret colour bimodality as a manifestation of two components rather 
than two distinct galaxy populations. Note that we have shown the 
original E/SO classifications without the additional colour tinkering 
described in Section 3.1 and the BS population shows up as fairly 
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Figure 10. {upper left) the observed {u — r)g versus Mb^^qq distribution (dots and 0.2 dex contours) and histogram {upper right), {lower left) the volume- 
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distinct and lying within the blue peak. Fig. 1161 (right) shows the 
observed distributions for the MGC continuum classifications of El 
(upper), Sa (middle) and Sc (lower). This suggests that while the 
El-Sa division is fundamental the Sa-Sc division is probably not, 
but instead merely an arbitrary cut of a continuous distribution. 



5 DISCUSSION 

Our main result is that bimodality occurs in both the rest-(u — r) 
and log(n) distributions and Fig. llSl highlights the 2D nature of 
this bimodality with two peaks of approximately equal height and 



shape in the observed (Af_B„(-p < —16 mag) plane (upper). We 
have argued in the previous section that these two peaks are best 
explained by the two component nature of galaxies, i.e., bulges and 
discs. Fig |17l shows how the stellar mass is distributed across the 
(u — r)c — log(n) plane today, and the red concentrated peak is 
noticeably larger containing (~ 71 ± 6) per cent of the total stellar 
mass. Spheriods (i.e., E/SO[red], c.f., Table|2} account for (35 ± 2) 
per cent. Systems with both bulges and discs must therefore ac- 
count for the remaining 36 per cent. If one naively assumes a mean 
B/T ~ 0.5 ± 0.2, then this yields a total stellar mass content in 
spheroids and bulges of (53 ± 7) per cent. The discs hence make 
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Figure 11. {upper left) the observed (u — r)c versus Mb^^qq distribution (dots and 0.2 dex contours) and histogram {upper right), {lower left) the volume- 
corrected {u — r)c versus Mg^^^ distribution (dots and 0.2 dex contours) and histogram( tower right). 



up the remaining (47 ± 7) per cent. Note that aUhough we truncate 
our distributions at Mbugc ^ 5 log ft = — 16 mag, the luminosity 
functions are sufficiently flat to suggest a minimal contribution to 
the stellar mass density and therefore dwarf systems are ne gligible 
in ter ms of total stellar mass estimates (as also inferred in lDriveil 
Il999h . 

The very narrow distribution of the red peak (see FigsQand|2j 
also implies a universal continuum shape which may only come 
about for old systems which formed with a universal IMF. Mergers 
which lead to bulges must either occur early or without significant 
star-formation i.e., dry-merging (see iGrahamI Eool Ipaber et all 



l2005l and lBell et alj2005l for fairly stringent constraints on the level 
of dry-mergers at z < 1). Hence it seems that around half of the 
stellar mass fraction (~ 53 per cent) was formed early with the re- 
mainder (~ 47 per cent) forming more recently. This at face value 
ties in with the well established decline in the local star-formation 
rate since z ~ 1 (Hookins 2004, Bell et al. 2005). However it is 
difficult to envisage intuitively how a 2D-bimodality could occur if 
the whole story consisted of a single and gradually declining star- 
formation rate. One can infer an approximate age for each peak 
from the MGC continuum fitting process in which the correspond- 
ing synthetic spectra ages were 15 Gyr(El), 7.4 Gyr(Sa) and 2.2 
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Figure 12. (upper left) the observed log(n) versus Mb^^qq distribution (dots and 0.2 dex contours) and histogram (upper right), (lower left) the volume- 
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(see Graham & Guzman 2003). 



Gyr(Sc). These equate to z > 3, z « 1 and z « 0.3. That the 
Sa and Sc systems both he in the same blue peak points towards a 
common and continuous formation process stretching and dechn- 
ing from at least z ~ 1 (i .e., the evolving star-forming population 
sought bv lBell et alJ200.^ . 

The red peak remains less certain, but because of the bimodal- 
ity would appear to share a distinct formation mechanism occur- 
ring at an earlier time (based on the narrower red colour distri- 
bution). That the core colour exhibits a narrower peak than the 
global colour suggests that it is the cores which are coeval rather 
than the entire red galaxy population. Of course red bulges can 



come about through early formation of sub-units and later dry- 
merging (i.e., a merger event in which minimal star-formation 
occurs due to the already depleted gas reservoirs). Recent con- 
straints on the levels of dry-mer ging appear th i s process occurs 
at a low level at z < 1 (see iGrahanj |2004 iBell et alJ Eoo3 
and also iFaber et aljEoOSh . Perhaps more persuasive is the ar- 
gument that many dry-mergers would eliminate the colour and 
metallicity gradients that one sees in early-type systems (e.g., 
iLa Barbera et alj|2003l lOe Propris et aljEoOS) . Carrying forward 
the notion of a single coeval event with an implied formation peak 
at z > 3, the most obvious connection to make is with the peak 
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Figure 13. {upper left) the observed po versus Mb^hqi- distribution (dots and 0.2 dex contours) and histogram (upper right), (lower left) the volume-corrected 
Po versus Mg^^^ distribution (dots and 0.2 dex contours) and histogram (lower right). The sohd line shows the Imown relation for spheroids (see Graham & 
Guzman 2003). 



in the quasar luminosity density at 2 « 3 fran et al."2004|). The 
MsM BH — L,a,n relations ( Novak. Faber & Dekel 2005) pr ovide 
additional support for this connection (see also Silk & Rees 1998 
for a more theoretical basis for this connection). The key question 
is the n whether the bulges formed through a s ingle monolithic col- 
lapse iEggen. Lyndon-Bell & S andage"l962h process, or through 
the rapid me rger of high-m ass components in high-density environ- 
ments (e.g.. lMenci e t al. 2005). An interesting observational con- 
straint comes from the phenomenon of 'core-depletion' in which 
some small fraction (~ 0.1 per cent) of the stellar population 
in the core region may be missing due to ejection during SMBH 



coalescence. Recent constraints jGrahadliool imply that giant 
spheroids (Mr ~ —22.5 ± 0.5 mag), have typically undergone 
one dissipationless major merger (involving the coalescence of two 
SMBHs) . If universal, this may point towards a n early monolithic 
collapse iEggen. Lvndon-Bell & Sandagell963) with no more than 
one major merger. Certainly any merging scenario would also 
need explain how th e Afgjyjgjj — a iF errarese & Merritt 20o3 
iGebhardt et alj2000l) and i\/sMBH ~ " iOraham et al. J200lh re- 
lations (see Novak et al. 2005) are preserved. However simulations 
involving multiple mergers occurring over the full lifetime do ap- 
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Figure 14. (upper left) the observed HLR versus Mb^^qq distribution (dots and 0.2 dex contours) and histogram (upper right), (lower left) the volume-corrected 
HLR versus Mg^^^^ distribution (dots and 0.2 dex contours) and histogram (lower right). 



pear t o be able to produce both the col our bimodality iMenci et alJ 
120051) and preserve the above relations iCroton et allZOOsI) . 

Based purely on the empirical data alone, the above discussion 
suggests that the two peaks are due to bulges and discs (rather than 
two galaxy populations) which may be associated with the peak in 
quasar activity and the broader star-formation distribution respec- 
tively which exhibit distinct trends over distinct time-scales. First 
we have the formation of the AGN/SMBH/bulge trinity, peaking 
at 3 followed by a more extended second phase of disc for- 
mation (presumably through splashback, infall and further accre- 
tion). To pursue this avenue further and in a rigorous rather than 



speculative manner now requires accurate bulge-disc decomposi- 
tions. These are currently in progress and results will be reported 
in forthcoming papers by Allen et al. and Liske et al. 



6 SUMMARY 

We have used the Millennium Galaxy Catalogue to explore vari- 
ous traditional methods for subdividing the galaxy population. The 
parameters include eyeball morphological type, fitted global Sersic 
index {n), 2dFGRS ?7-parameter, rest-(it — r) colour, MGC con- 
tinuum shape, extrapolated central surface brightness, half-light ra- 
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Figure 15. (upper) the observed distribution of log(n) versus (u — r)^ colour (left, shown as data points and 0.2 dex contours; right shown as a 3D plot with 
the z-axis linear), (lower) the volume connected distribution of log(n) vs {u — r)c colour (left, shown as 0.2 dex contours and right as a 3D plot with z-axis 
linear). 



dius and stellar mass-to-light ratio. In both rest- (it — r) colour and 
log(n) the observed and volume-corrected distributions are clearly 
bimodal and the joint distribution show obvious 2D-bimodality. We 
find that morphological classification is susceptible to contamina- 
tion from smooth, blue, low-luminosity systems and that overall 
the traditional morphological system is wholly problematic to im- 
plement. We find that continuum fitting of the broad-band colours 
leads to a good separation of the two populations which we label 
"old" and "young" (and corresponds to the "re d" and "blue" peaks 
respectively, identified by Strateva e t all200ll) . 

The E/SO systems lie predominately in the old, red peak and 



the Sd/Irrs lie wholly in the young, blue peak with the Sabc class 
straddling both peaks and the divide. We take this as strong evi- 
dence that the galaxy population does not consist of two classes 
but two components, consistent with the classical idea of bulges and 
discs. The red peak hence constitutes older bulges and the blue peak 
younger discs. We advocate that this joint-bimodality (i.e., colour 
and concentration), may reflect two distinct formation mechanisms 
occurring over different time-scales. We suggest that the red peak 
may be associated with AGN/SMBH formation via monolithic col- 
lapse and that discs come about through later secondary processes 
(splashback, infall and accretion). Forward progress on galaxy for- 
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Figure 16. (left side) The distribution of E/SO (upper), Sabc (centre) and S 
(upper), Sa (centre) and Sc (lower) in the log(n) vs {u — r)c colour 

mation demands routine bulge-disc decomposition to enable the 
study of these potentially distinctly formed components. This is 
currently in progress. 




■s (lower) in the log(n) vs {u — r)c colour (right side) The distribution of el 



ucts are publicly available from |http://www.eso.org/'^jliske/mgc7| 
or on request from J. Liske or S.P. Driver. 



The Millennium Galaxy Catalogue consists of imaging data 
from the Isaac Newton Telescope and spectroscopic data from the 
Anglo Australian Telescope, the ANU 2.3m, the ESO New Tech- 
nology Telescope, the Telescopio Nazionale Galileo, and the Gem- 
ini Telescope. The survey has been supported through grants from 
the Particle Physics and Astronomy Research Council (UK) and 
the Australian Research Council (AUS). The data and data prod- 
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Table 2. The Millennium Galaxy Catalogue S-band luminosity functions for the overall sample (all) and various subdivisions (as indicated). Errors are purely 
Poissonian. 
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